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Abstract 

We first study dark energy models with a minimally-coupled scalar field and exponential poten- 
tials, admitting exact solutions for the cosmological equations: actually, it turns out that for this 
class of potentials the Einstein field equation exhibit alternative Lagrangians, and are completely 
integrable and separable (i.e. it is possible to integrate the system analytically, at least by quadra- 
tures). We analyze such solutions, especially discussing when they are compatible with a late time 
quintessential expansion of the universe. As a further issue, we discuss how such quintessential 
scalar fields can be connected to the inflationary phase, building up, for this class of potentials, a 
quintessential inflationary scenario: actually, it turns out that the transition from inflation toward 
late-time exponential quintessential tail admits a kination period, which is an indispensable in- 
gredient of this kind of theoretical models. All such considerations have also been done including 
radiation into the model. 
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I. INTRODUCTION 



Scalar fields have been used extensively in cosmology since inflationary theories were first 
conceived in 1981 [l|. They are in fact useful to fill in the universe with a negative-pressure 
content, hence giving the typical inflationary behavior, with its exponentially increasing scale 
factor a{t) ~ exp{HQt), where Hq is the present value of the Hubble parameter. Scalar fields 
have been used both in minimal and non minimal coupling with geometry. (An overview of 
scalar-tensor theories is given in Refs. ; a more updated review can be found in Ref. 

When observational results on Type la supernovae (SNela) have stronglysupported the 
possibility that the universe is now in an accelerated stage of expansion js] jg] j3] jsl, it 
has become necessary to consider again a negative-pressure component. The late-time 
infiation so generated is however soft with respect to the earlier one and is again based on 
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the dominance of at least one scalar field referred to as quintessence 19 1 Il0| 

n n 

or, earlier, as x-field [jj [15l. This gives rise to a sort of energy component which only 
appears measurable through gravitational effects. There is in fact a wide agreement that 
today a form of such a dark energy |l^ 1^1 has to be taken into account in any realistic 
— mode,. Seen . a constant A-tpn Q fl, dark enetgy den.ty has been 
usually seen as the vacuum energy density [21|, leadin g th erefore to relevant problems due 



to huge discrepancies between theory and experiments 22|] 19|] 20|] 23 1. 



In this context, the scalar field in itself has gained new attention, and speculations about 
its nature and evolution have got new strength. Papers like Refs. 2J] and 25|| have been 
reconsidered again under a new light, as they show much of the role a scalar field could 
play when one wants to have an accelerated expansion today. Many models have thereof 
been constructed, first of all trying to consider more appealing kinds of potentials driving 
the dynamics of the scalar field. Among others, the potential has in fact to be seen as 
an important ingredient of all the involved theoretical models, and many are the trials 
to give prescriptions in order to reconstruct its form according to observational data (see 
's. |26 [_ 27|, f or instance; but, for more theoretical considerations, see also Refs. 28| 



12| [29| [30j). Here we take into account a specific class of potentials useful for 
the present evolution of the universe, later investigating all the relative involved cases so 
derived in the cosmological equations. The exponential-type potential has received much 
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attention since the late 80's 



3l| |32| y , (but see also Refs. 



33[ [3^ [35 
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39[ |29| [40|, [4l|], just to cite only some other old papers). In Refs. [4^J [4a|J we 
have found general exact solutions for two classes of exponential potentials for a scalar 
field minimally coupled to gravity, in presence of another dust component (ordinary non- 
relativistic pressureless baryonic matter), also showing 4^ that such solutions seem to fit 
SNela data given in Ref. Gj. (See Ref. 45| for much more on confrontation of theoretical 
predictions with observational data.) Furthermore, in Ref. 46|] other approximate and exact 
solutions for cosmology with exponential potentials can be found, while up to more recently 
it is still widely referred to (as, for example, in Ref. 4^, where the Noether Symmetry 
Approach ^ is used to probe the nature of dark energy, well underlining the role of 
the exponential potential in such a task). This seems to us a good motivation to investigate 
also on generalizations of this kind of potentials. (As a matter of fact, it is often studied 
also the case with the hyperbolic sine potential, which can easily be assimilated with some 
situations studied below). 

Moreover, also for such a class of generalized exponential potentials it is possible to 
exhibit some special properties: it turns out, actually, that the Einstein field equations are 
integrable (in the Liouville sense) and separable, that is they can be analytically integrated, 
at least by quadrature. Finally, it is worth noting that our class of exponential potentials 
can be selected by finding the most general variables transformation which diagonalizes the 
scalar field kinetic-energy form, leaving the transformed Lagrangian simple. 

Therefore, among other considerations, we then generalize the exp onential potential to 
the extent allowed by the particular technique used in Refs. 4^ 431] in order to integrate 
cosmological equations, still bearing in mind that, in order to describe the present universe, 
solutions should always allow a late-time acceleration. Such a technique is generally known 
as the Noether Symmetry Approach to cosmology j48|, j^l. (Note that in Refs. {3], js], such 
an application of the Noether Symmetry Approach already leads to prefer the exponential- 
type potentials.) As in Refs. j42] and j43| we anyway do not completely adopt it here, 
limiting ourselves to import from it only the change of variables needed to solve equations. 
As a matter of fact, we have to point out that this transformation (together with the choice 
of a specific kind of potential) has to be considered as the main drawback of the approach 
above, since it leads to a point symmetry which probably applies only to the background 
evolution, having nothing to do with the actual symmetries of the full theory. We analyze 
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the exact solutions we find in this way, discussing when they are compatible with a late time 
quintessential expansion of the universe. 

As a further step, we also consider and illustrate a possible scenario in the framework 
of the quintessential inflation paradigm with scalar fields, where an inflationary potential 
drives also the quintessential phase of the scalar field evolution, by means of an exponential 
form of its late time tail. In this connection, we discuss how such an evolution mechanism for 
the scalar field potential can be compatible with the so-called inflation-kination transition, 
when the field energy density is dominated by the kinetic energy of the scalar field (f. All 
this clearly leaves apart the greater complexity involved by the more general forms of the 
potential studied above, and will be the subject of a forthcoming paper trying to discuss 
more in general the quintessential inflation. Even if in this second paper we should be able 
to overcome the dichotomy present in this work, where the two different parts of the paper 
could well live apart, we believe that it is nonetheless already interesting to offer both the 
issues in a single paper, trying to unify the whole subject a little forcibly. 

In Sec. II we introduce the class of potentials and the general cosmological setting for 
further considerations. In Sec. Ill, we systematically derive, when possible, general exact 
solutions case by case. In Sec. IV we discuss connections between exponential potentials 
quintessence and inflation. Finally, Sec. V is devoted to a conclusive discussion. 



II. COSMOLOGICAL MODELS WITH (GENERALIZED) EXPONENTIAL PO- 
TENTIALS 



Motivated by the considerations above, let us continue to assume a Friedmann-Lemaitre- 
Robert son- Walker (FLRW) metric and fix the curvature scalar k = 0, since the spatially 
flat situation appears to be the most appropriate one according to the CMBR observational 
data (see Refs. 49| 50| for recent information). In what follows we have to consider the 



period of life of the universe after the decoupling time, in order to be realistic when taking 
the two components like dust and scalar field ip into account, like in Refs. 42| 43 1. 

If if is minimally coupled to gravity, and considering no special choice for constants, the 
cosmological equations are written 

SttG 



3H' 



-ipm + Pv) 



(2.1) 
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H + H'' = -^{p^ + p^ + 3ipm + P^)), (2.2) 

>^ + 3Hip + V'{ip) = 0. (2.3) 
The fluids filling the universe have equations of state given by 

Pm = , p^ = w^p^ , (2.4) 

being pm = Da~^, where the parameter D = pmoOo^ is determined by the current values of 
Pm and a. Let us also recall that 

= + V{^) , ^ - V{v) , (2.5) 



and 



_ - 2V(,p) 



The cosmological equations can then be rewritten as 



SttG , 1 



a J \ 2 



a f SnG ( 1 



Da' + -0^ + V{ip) , (2.7) 



nn 

It can be shown |48|],[3[ that the last two equations can also be deduced from an action 
principle based on the point Lagrangian 



(2.10) 



In this way, cosmological dynamics can be considered on a space with two coordinates a and 
(fi, a and (p being the velocities j^sj j^. The fact that C has a constant additive term is 
understood by considering Eq. (12. 7p . which can be seen as Ec = 0, where Ec = |if + ~ 
is the so called energy function associated with C. We find the constant to be zero for 
physical reasons (i.e., induced by the homogeneous and isotropic limit of Einstein's general 
field equations). In this way, when Eqs. (12. 8 p and (12. 9p are solved, Eq. (12. 7p is nothing but 
a constraint on the integration constants involved. 
On defining 



T = 3aa'-^a'^\ U = -^{a'V{v) + D) , (2.11) 
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Eq. fl2.10p formally becomes C = T — U. This has been already noted in Ref . 5l| , where 
we found the most general transformation 



a = f{z, w) ip = g{z, w) (2.12) 
leading to a diagonalized and simpler form of the kinetic energy 

r = a^z^ - (3^w\ (2.13) 

with a and (3 nonnegative real numbers. This simplification in T implies, of course, a 



complication in the transformed U'. However in Ref. [48|] we also found that there is at 
least one class of potentials not giving such a complicated expression for U' , rendering, 
on the contrary, the transformed C simpler than C. It turns out that such potentials are 
exponentials: 

V{ip) = Vo {A^ exp (2C^) + exp (-2C^) - 2AB) , (2.14) 
with C = (3ttG)/c'^ and A, B real parameters. They in fact allow an exact integration of 

nn 

the cosmological equations |48[ p[. Moreover, it is worth noting that, for such a class of 
exponential potentials, 

y (y,) = Vi exp {2Cip) + ^ exp (-2C(^) + Aq , (2.15) 

being Vi, V2, Aq all free parameters, the Einstein field equations admit alternative La- 
grangians, which is a circumstance somehow exceptional for a dynamical system, with many 
meaningful consequences: actually, it turns out that it is integrable (in the Liouville sense) 
and it is separable, i.e. there is a suitable change of variables by which it is splitted into 



separated one-dimensional systems (which are integrable by quadratures j52j). 

We thus feel motivated to continue considering exponential-like potentials, taking now 
the class of potentials 

V{<^) = A^ exp {a<^) + eB'^ exp {-a<^) + A , (2.16) 

where e = ±1, a = ^/VlixG j (? is a fixed constant. A? and B^ are arbitrary nonzero param- 
eters, and A ^ 0. Eq. fl2.16p generalizes the exponential potentials already considered in 



Ref. 42|, which are 

V[lp) = B'^fty.^{-aLp) , 1/((^) = exp ((T(^) + exp (-a(^) . (2.17) 



These cases (where e = +1 and A = 0) are therefore omitted in this paper. Moreover, as 
noted, the case with e = +1 and A = —2AB has been already treated and discussed in Refs. 
Ql ^1 S) cind it will not be touched upon again here. 

Since the theory is invariant under if — > — the class of potentials in Eq. fl2.16p is 
the most general one of this type and also includes the first one shown in Eq. fl2.17p . 
V{ip) = B"^ exp {—(Tip) is in fact derived from the expression in Eq. f l2.16p just setting 
(besides e = +1 and A = 0) also = 0, or 5^ = and ip — )■ —ip. (Of course, in order to 
get exact general solutions, the treatment is different from the beginning for each one of the 
potentials in Eqs. fl^T^ and I KTT^ .) 

As (among others) in Ref . 42| , then, let us use again the transformation 



n — 1 , 

a = , = - log 

4 a 



Biu + v) 



(2.18) 



A{u - v) 

leading to the new variables u and v. Such a change of variables is invertible, provided that 
a 7^ 0. Of course, when the solution for the scale factor eventually has a zero in the future, 
we cannot follow it beyond that point. We must in fact consider such a solution as good 
just between this zero and the eventual other one it can have in the past. (Later on, we 
will comment on this again.) Since we also want a > 0, Eq. f l2.18p imposes the obvious 
restriction > v"^, or {u + v){u — v) > 0, giving u > v when both u and v are nonnegative. 
This soon implies that A and B must have the same signs (i.e., AB > 0). In what follows, 
we choose to consider A > 0, B > 0,u > v, and m > 0, f > 0. We also disregard the absolute 
value sign in the expression for ip, which does not really limit the class of potentials we are 
studying. By virtue of Eq. fl2.18p . the potential in Eq. fl2.16p becomes 

+ (2-19) 

u — V 

and we define 

V+{u, v) = 2AB——- + A (2.20) 
u — V 

for e = +1, and 

uv 

V4u,v)=AAB- - + A (2.21) 

— t)^ 

for e = — 1. 

In the case when A'^ = 0, anyway, we have to use the change of variables 

1 u 

= uv, ip = log - , (2.22) 

a V 

7 



as already done in Ref. 42| when A = 0. This case will be treated in detail in Sec. Ill, 
taking into account the possible values of A. 

Now, transforming the other terms in Eq. f l2.10p by means of Eq. fl2.18p , we can eventually 
arrive at two expressions for the point Lagrangian, according again to the two opposite values 
of e, 

2 

C+ = u^ -v"^ + ^ [{2AB + \)u^ + {2AB - \)v'^ + AD] , (2.23) 

2 

C_=u^-v'^ + ^ [AABuv + \{u^ - v"^) + AD] (2.24) 

where the possibility of multiplying the Lagrangian by a costant, equal to 3 in our case, has 
been exploited In order to find equations and solutions for the cosmological model, at this 
point, it is better to consider the various different cases separately, and this is what we are 
going to do in the next Section. 

In the following, we fix four conditions (see Ref. 43|] for further details). First of all, we 
set the origin of time by choosing a(0) = 0, paying attention to the fact that, even if this can 
be done generally, it is not a real initial condition for a{t). The initial time tin chosen for 
any specific model is in fact arbitrary, to some extent, and mathematical solutions can also 
present more than one zero, eventually infinitely many. In such cases, one must be careful 
in choosing only one physically meaningful period of time for the evolution of the universe, 
noting that we can accept no more than one zero at the beginning of time. An eventual 
second zero occurring in the future has then to be recognized as a time of recoUapse of the 
universe. What happens later thus loses any meaning. We have thus to set m(0) = or 
f (0) = in evaluating initially Eq. fl2.22p. b eing it arbitrary for the moment to decide which 

42| . for instance, we decided to set to zero both, so 



one is better on solid ground. In Ref. 
selecting the tracking solution. 

The second condition that we assume is that the present time (age of the universe) is the 
unit of time, to = 1- Being tin unknown, this is not exactly the age of the universe, but the 
difference can be considered irrelevant for our purposes; anyway, even if it were be possible 
to avoid this condition, we thus get rid of a badly known quantity. Our third condition is 
to set Oq = a.(to) = c^(l) = !> which fixes the normalization of a as standard, while the last 
condition is to set H{tQ = 1) = 'Hq. After the choice of to? this latter parameter turns out 
to be of order 1, even if it is not the same as the usual h. As a matter of fact, since we 
are using an arbitrary unit of time, such a parameter does not give any information on the 
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observed value for Hn 



43|. 



As to A = 0, as already said, some considerations are given in Ref. 



42|. When A = = 



and e = +1, for instance, Eq. fl2.22p easily leads to general exact solutions for a(t) and ip(t) 



42]. For A = -B^ = 0, everything works the same, since there is symmetry in the potential 



in Eq. f l2.16p with respect to a change of sign in (p. The situation where 7^ 0, -B^ 7^ 0, 
and e = +1 has also been partially examined in Ref. |42|. 

In the next Section, we will thus focus our attention mainly on the case with 7^ 0, 
B"^ ^ and e = — 1, choosing not to consider at all the situation with X = = and 
e = — 1, since it involves an always negative potential and may also give a negative energy 
density, by virtue of Eq. f l2.5p . Most of all, let us note that it never allows the pressure to 
be negative, which forbids the possibility of describing the accelerated expansion we observe 
today. 

We will also rule out the trivial case given by = A, while the remaining situations will 
be dealt with systematically. 



III. SOLUTIONS 

Once we have fixed our Lagrangian expressions in Eqs. f l2.23p and fl2.24p . we can soon 
derive the related equations for cosmology as Euler-Lagrange equations. If our transforma- 
tion in Eq. f l2.18p then works, such equations should turn out to be solvable, even if, of 
course, exact integration is not always easy. The analysis strongly depends on the relative 
values assumed by the constants involved. Thus, in what follows we separately discuss each 
situation generated, as a first step, by a different choice of A values. Then, we pass to further 
investigate the features of the equations and their solutions, starting from consideration of 
the two allowed values of e. 



A. The A = case 



As mentioned, this situation has already been partially treated in Ref. [42j, where we 
consider only three cases: i) = and e = +1, ii) = 0, and iii) A^ ^ 0, B"^ ^ 0, 
and e = +1. The case with A^ = and e = — 1 has already been touched upon in Sec. 
II, noting that it is not so interesting here in our considerations on present cosmology. On 
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the other hand, when both and B'^ vanish, the potential vanishes and (f is free, giving 
= = 0^/2 > 0, and w^p = 1. This means that the scalar field behaves like stiff matter 
and introduces a term oc in the cosmological equations, thus becoming soon negligible 
in the expanding evolution of the universe. In this way, ip would never produce a late-time 
inflation, for instance, as instead recent observations seem to require. 
There only remains one case to discuss, i.e.. 



-1 



(3.1) 



With such assumptions, setting A = in Eq. fl2.24p yields 

2AttG 



-{ABuv + D) 



(3.2) 



so that the related Euler-Lagrange equations are 



u = uj'^v , V = —uP'u 



(3.3) 



where 



VI-rGAB 



2_ ..n^^^ =^-^AB 



It can be shown that the solution reads as 

uit) -- 



C3 + C4 1 / X , . .i^t X 

cosh(^p) + C2 smh(^pj 



W2' 



Lot 



sm 



V2 \ W2 



C3-C4 . ,,ujt 

— — — smhl^^i + Ci coshl^^ 
. V2u W2 V2 



cos 



. ut 

71 



(3.4) 



vit) 



C3 - C4 , , . ut 

■ cosh(^p) — ci smh(- 



V2i 



W2' 



ut 



sm 



72 J V2 



+ 



C3 + C4 . ,,ujt ut 



ut 



cos(^) , (3.6) 



so that we get two functions of time t depending on four arbitrary parameters (mi, M2, ti, 
and ^2)- By virtue of the energy constraint Ec_ = 0, which is written 

247rG' 



• 2 -2 o 2 

U — V — 2u UV — 



-D = 0, 



(3.7) 



Eq. (1331) yields 



D = AB{^-^-c^c, 



(3.8) 
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If we define = (cs + c^)/ {\/2ll!) and = (c3 — Ci)/{\/2uj) we then get D = 745(^3^4 — 0102). 
The initial scale factor at t = is such that the volume a^{t) is at the beginning a^(0) = 
(ci^ — C2^)/4, in light of Eq. (I2.18p . For sake of simplicity, we can fix the origin of time in 



such a way as to get a vanishing scale factor (as in Ref. 42|]). This, of course, is completely 
arbitrary, and we must always remember that the real beginning of time for the model 
under examination is actually afterwards, i.e. , at an after-decoupling moment which can be 
arbitrarily delayed for the kind of considerations we are making here. This also means, then, 
that we are setting ci = C2 or, possibly, ci = C2 = (since we have put a(0) = 0). Anyway, 
both choices still allow D to be nonvanishing and positive, which we certainly require. Of 
course, this would be gained even taking C2 = simply. But let us here put Ci = C2 = 0, so 
that we can write 

D = ABk^ki , (3.9) 

from which we find that the two constants k^ and k^ have the same sign, i.e., k^k^ > 0, 
implying that Cs^ > c^^. 

Let us also note that, from Eqs. (13. 5p and (13. 6p . we get u(0) = — ci and t>(0) = C2. Thus, 
independently of all considerations above, Eq. (12.18^ 9 implies 

1 BiuiO) + vm 1 B{-c, + C2) 

= - log -— = -log— -, 3.10 

a A[u[0) — v{0)) a A(— ci — C2) 

and the choice Ci = C2 or, equivalently, a(0) = 0, soon gives ipiOi) = 0. 

By substituting Eqs. (13. 5p and (13. 6p into the expressions in Eq. (I2.18P for a^(t) and f(t), 
after some algebra we find the cosmological solutions 



a{t) 



Q (/cs^ - k4^) (^1 - cos{V2ujt) cosh{V2ujt) 



+ 2k3k4sm{V2ujt)smh{V2ujt)\y^^ , (3.11) 



l/a 



^ [b [{k^ - ki) sin(y2wt) cosh{^/2ut) + {k^, + k^) cos{V2ut) sinh{^/2ujt)] \ 
^ ^^\A [{k^ + ki) sm{y/2cot) cosh(v/2a;t) + {k^ - k^) cos{y/2ujt) smh{y/2ut)] J 

(3.12) 

The scale factor in Eq. (13. lip yields a sort of cyclic universe with spatially flat hyper-surfaces 
(being k = 0). Since, in practice, it is nonvanishing only in an infinite number of limited 
intervals, this means that the solution describing expansion should be considered only in 
one stage (in its first stage, for example). 
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On the other hand, the present value of the matter density parameter is given by 

n^o = 8nGD/{ao^Ho'^) = StiGD/Uq^ ~ SttGD = SuGpmo , (3.13) 
being = 1 and Hq = ~ 1. 

B. The A > case 

Here, first of all, we have to distinguish between the two cases with e = +1 or e = — 1, 
i.e., between Lagrangian functions in Eq. fl2.23p or Eq. fl2.24p . respectively. Then, as we 
can see, our considerations must also take the relative values of AB and A into account. 
With respect to this, note that the situation given hj A"^ = B"^ = describes a universe 
filled in with dust and free scalar field in presence of a positive cosmological constant, since 
the potential V{ip) reduces to A only. This means we shall limit our considerations here to 
AB ^ 0. 



1. The € = +1 value 



First of all, let us not e ag ain that the case with e = +1 and A = —2AB has been 
already studied elsewhere 48] jsil ^. Thus, we do not comment on it here. Now, on the 



a 



other hand, bearing in mind that u"^ > v"^ and AB > 0, the choice e = +1 always gives a 
positive potential for the scalar field and leads to two different cases, since we can write (for 
2AB ^ A) 

(3.14) 
(3.15) 

(3.16) 
(3.17) 



-{2AB + X) =ui^ 
{2AB - A) = ±0/2^ 



so that Eq. fl2.23p becomes 



We thus deduce the equations 

il = OJi'U , V = ±U]2^V , 



where the plus (minus) sign corresponds to the minus (plus) sign in the Lagrangian In 
both cases we find 

u{t) = a exp (cjit) + (3 exp {—uit) , (3.18) 
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a and /3 being integration constants. Furthermore, we find 

vii) = Vx^viiiuj^t \ V'l) (3.19) 

when the plus sign {2AB > A) is taken in the Lagrangian, or 

v{t) = f 1 exp {uj2t) + V2 exp {—uj2t) (3.20) 

when the minus sign {2AB < A) is chosen there, {vi and V2 are integration constants.) 
The integration constants are constrained by means of the equation Ec_^ = 0, which is 

■ii'^ -v^ - ui^u^±U2^v^ - 2a^D = . (3.21) 

This gives a relationship that makes it possible to determine the more physically meaningful 
constant D in terms of the other constants, i.e., 

D = _(l2^t^±±^ > (3.22) 

when the minus sign is taken in Eq. (13.211) . and 

D = : "^"^'^ > (3.23) 

with the plus sign. Focusing on the situation described by Eq. f l3.19p . i.e. choosing 2AB > A, 
we see that the inequality (13.221) implies that — 4a/3a;i^ > t>i^C(;2^, that is, af3 < —vi'{2AB — 
\)/['i{2AB + A)], giving a/3 < 0. So, when (3 = —a, one finds u(t) = 2a;sinh (oJit), and we 
may write the volume a^{t) as the difference between the squares of a hyperbolic sine and a 
sine, where now Aa'^ui^ > Vi'uj2^- On the other hand, choosing f 2 = does not change the 
sign of the value of D and implies a(0) = 0, i.e., we can fix the origin of time as before, as 
well as in Ref. [42]. Of course, putting w 1 = also leads to a(0) = 0, but then we always 
have v{t) = 0, so that a^{t) = sinh^(a;it), being m^(0) = 0. Anyway, let us stress that, in 
the case under examination, it must be a(0) = in order to avoid an initial negative value 
for the scale factor. 

At the same initial time t = 0, the scalar field generally has the constant value 

1 (3.24) 

a A[a + p — vi sm f 2) 

which gives rise to an undetermined form for f3 = —a and V2 = 0, unlike what was found in 
the case studied above, when A = 0. 
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The cosmological solutions can be then written as 

"4a^ sinh^ {uit) - vi^ sin^ (cu2t) 



a{t) 



4 



(3.25) 



= i log ^[2asinh(a;,t) + .,sin (c.,t)] 
^ (T ^ A[2a sinh (wit) - sin (c^at)] 
so that the asymptotic behavior of the scale factor is an exponential one. As before, it 

thus expresses a late-time nonsoft accelerated expansion of the universe, but a better agree- 
ment with observations demands other considerations as to a more refined comparison with 
astrophysical data. 

On the other hand, discussing the other allowed situation, with 2AB < A, the condition 
f l3.23p implies that viV2UJ2'^ > that is, a/3 < ViV2{X — 2AB) / {2AB + A), and now we 

can only say that f if 2 < gives a/3 < 0, while f if 2 > leaves the possibility of having also 
a/3 > 0. But let us note that a{t) > involves m — f > and, therefore, m 7^ 0, which means 
that a(3 7^ 0. If we then put (3 = —a (i.e., a/3 < 0), we get u{t) = 2asmh.uit, as before; 
this, in turn, requires that a^ > — fif2(A — 2AB) / {2AB + A), giving no restriction on the 
sign of V1V2, of course. 

Thus, we can also choose V2 = — fi (implying a^ > fi^(A — 2AB) / {2AB + A)) and get 
v{t) = 2vi sinh {u2t) , so that we may write the volume a^{t) as the difference between the 
squares of two hyperbolic sines. As before, this choice implies a(0) = for every value of a 
and Vi. Moreover, at t = 0, the scalar field has in general the constant value 

ip{0) = - log -, 3.27 

a A[a + p — Vi — V2) 

which is again undetermined when we choose /3 = —a and f2 = —vi- 
Now, the cosmological solutions are written as 

a{t) = [a^ sinh^ (uit) - vi^ sinh^ (^2^)] , (3.28) 
1 B[a sinh (uit) + f 1 sinh (^2^)] , . 

V2(t) = — log— 7 — — r — — TY . (3.29) 

cr A[a smh [uitj — Vi smh [uj2t)\ 
Of course, the expression of a(t) is non negative only when the constants a, ui, vi, and U2 

are such that the difference present in a{t) is positive. In such a case, when U2 < u)i we have 

an asymptotic exponential regime for the scale factor, as the one we already found before, 

with the other choice for the signs of the constants. The case with 002 > Wi, on the other 

hand, has to be discarded since it never happens. It asymptotically leads to a non physical 

situation, with a negative infinite value of the scale factor. 
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2. The € = —1 value 



Let us consider, now, the Lagrangian in Eq. f l2.24p . resulting from the potential V- 
expressed in Eq. fl2.2ip . Since has to be positive, as we have chosen, it must be 

AABuv > -A(u2 - v^) , (3.30) 

being u"^ > f ^, hence allowing uv ^ 0. (Note, in fact, that when uv = 0, the potential simply 
reduces to V- = \, therefore becoming negative.) 
Putting 

= (T^AB , = ^\a^ (3.31) 

into Eq. f l2.24p gives 

£_ = {{2 - {;2 + 2ui^uv + U2^{u^ - v"^) + 2a'^D , (3.32) 

so that the resulting Euler-Lagrange equations are 

u = oji'v + (jj^^u , V = —uji'u + uj2^v . (3.33) 

These equations constitute a linear system of homogeneous coupled ordinary differential 
equations of second order, which can, more appropriately, be rewritten as: 



(3.34) 



u{t) \ _ i dii ai2 \ I u{t) 
v{t) j y aai a22 j y v{t) 

where an = 022 = 1^25 ^^^1 012 = ul = —021. The method of solving such a system reduces, 
in fact, to diagonalizing the matrix 

A^h'^A. (3.35) 

\a2i a22 j 

Actually, if we introduce the rotation matrix R (whose column are the eigenvectors ipi,ip2 
of A), the transformation 

y = Rx (3.36) 

(where we have posed y = {yi(t) , y2{t)) and x = {u(t),v{t))) decouples our starting system. 
One has 

^=MAM~V, (3.37) 
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(3.42) 



where M.AR ^ is the eigenvalues' diagonal matrix. Such eigenvalues are 

-El = i (^-a/ (an - 022)^ + 4ai2a2i + an + ^22) , (3.38) 

^2 = ^ (^/(aii - 022)^ + 4ai2a2i + an + a22) , (3.39) 

and the eigenvectors associated with Ei and E2 are written as 

( a/ (an - a22)^ + 4ai2a2i - an + a22 ) 

^1 = ^ -e\ + 62 , (3.40) 

2a2i 

( a/ (an - a22)^ + 4ai2a2i + an - a22 ) 

^l'2 = ^ ^ei + 62 . (3.41) 

In terms of the transformed vector y, the system in Eq. fl3.33p reduces to a linear system of 

homogeneous decoupled ordinary differential equations of second order, which can now be 
integrated exactly: 

'yi{t) \ ^ ( El W yi{t) 
Ut) I \ ^2 y \y2it) 
The solution, in its general form, is 

l/i = cie^* + C2e-v^*, (3.43) 

1/2 = cae^* + C4e-^* . (3.44) 

The transformation Mr^y provides the analytical expressions of u{t) and v{t): 

u(t) = -U (^cieV-f-^-i + c2e-V-?-*-i + ie-V-?+^-i {c^ + cge^V-^+^-i j j , (3.45) 

v{t) = i (^ze-V^?-*^i (^C2 + cie^V^F^'^i) + cgeV^^+^'^i + 046"*^^^^) , (3.46) 

the Ci being (where i = 1....4) integration constants. It turns out that, in order to have a 
real- valued function representing the scale factor a(t), according to Eq. f l2.18p . we have to 
set C2 = C4 = 0, C3 = i, so as to obtain the following form of the scale factor: 

a{t) = ci^e2*V-t+-2cos(f) ^ ^3 47) 

where we can set Ci = juq, ag being a(to) and 7 being a positive constant. With such a 
choice of the integration constants we obtain a complex scalar field 



tsin (2t^wf + a;|sin(f)) cos f2t^Lof + Lol sin (|)) ^ 
0(t) = ^ ^ ^ ^ - - log (c?) + log {ool) . 

(3.48) 
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FIG. 1: Comparison between the time evolution of the real (solid lines) and the imaginary part 
(dashed lines) of the scalar field (j), for different values of the integration constants. It turns out that 
the imaginary part of the scalar field slightly oscillates around zero, and the real part dominates. 
The solid and the dashed lines correspond to C2 = C4 = , C3 = i , ci = 10 , /3 = arg + iuj'^ = ^ , 
uji = —7=- , and L02 = -7^- The solid red line and the dashed blue one correspond to n = 0.1, while 

\' 2n y In 

the others correspond to n = 0.05. 

Through an appropriate setting of the integration constants it is possible to see that the 
imaginary part of the scalar field slightly oscillates around zero, and the real part dominates, 
as shown in Fig. ([T]). In Eq. ( 13.47^ . on the other side, the scale factor exhibits an oscillation 
around an exponential growth for large t, hence offering the possibility to describe the now 
observed stage of dark energy dominance. 

C. The A < case 

As above, we have again to discuss separately the two cases with e = +1 and e = — 1, also 
considering the relative value of AB with respect to A. First of all, let us note that, when 
A? = B'^ = 0, we are once again in a universe with a free scalar field plus dust, but, now, in 
presence of a negative cosmological constant. Thus, we shall once again limit our following 
discussion to the situation with AB 7^ 0, still describing the analysis as systematically as we 
can, in order to offer a useful tool in the whole issue. 
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1. The e = +1 value 



Starting from Eq. (12.231) . let us remember that we are examining the case with AB > 0, 
that is, when a'^{2AB + A)/2 = iwi^ and a'^{2AB — A)/2 = U2^. This means that we have 
to consider two separate possible situations, due to the relative values of AB and A. 

The Lagrangian in Eq. f l2.23p is written as 

£+ = - v'^±uji^u^ + uji^v"^ + 2a'^D , (3.49) 
giving rise to the equations 

ii = iwi^M , V = —uj2^v , (3.50) 

where the plus (minus) sign corresponds to the minus (plus) in Eq. (I3.49p . The solution for 
u{t) is 

u{t) = uism.{(jjit + U2) (3.51) 
for the minus sign in Eq. fl3.49p (i.e., with 2AB < A), and 

u{t) = ui exp(ci;it) + U2 exp(— c^it) (3.52) 

for the plus sign (i.e., with 2AB > A), {ui and U2 are integration constants.) The solution 
v{t) is in both cases 

v{t) = asm{ijj2t + (3) , (3.53) 

a and /3 being integration constants. 

The constraint equation Ec_^ = now becomes 

- ± wi V - W2 V - 2(t2l) = , (3.54) 

leading to 

D = — > (3.55) 

for the plus sign in £+ (corresponding to the minus one in Eq. fl3.54p ). and to 

= - \^2 — > (3-56) 

for the alternative case with the minus sign in £+ (and the plus sign in Eq. f l3.54p ). 

From (I335I) . we must have U1U2 < -a^{2AB - A)/[4(2yli? + A)] < 0. Now, choosing 
U2 = —Ml, one finds u{t) = 2ui sinh(a;it), so that the volume a^{t) is written as the difference 
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(4Mi2sinh2(wit) - a'^sm\u2t + /3))/4, with u^^ > a'^{2AB - \)/[A{2AB + A)]. If we then 
set /3 = 0, this is consistent with the physically required positive sign of D and, at the same 
time, annihilates a(0). Note, on the other side, that the value a = gives v{t) = at any 
t, and a^{t) = ui^ sinh^(c(;it), so that this choice still yields a(0) = 0. If we do want to have 
both a(0) =0 and assume a = (leading to a finite constant <y9(0)as well as - as we already 
said-, as an initially vanishing scale factor), we instead (and unfortunately) find a constant 
scalar field at any time. Thus, even if it is true that, when both U2 = —Ui and /3 = 0, 
the argument of the logarithm in the <y9-field becomes an undetermined form, the condition 
a(0) = does require U2 = — mi and (3 = 0. 

With such choices the cosmological scale factor and the scalar field are 

a(t) = ('"''"""''"''^-°'""''"^") , (3.57) 

a A{2ui smh(u;it) — a sm[UJ2t)) 
But, even if this solution would be asymptotically good for describing an inflationary stage 
at present in the universe, we must rule it out because of the above considerations.. 

Turning, thus, to the situation leading to Eq. fl3.56p . we see that it has to be 
tii^ < a^((X'2^/u;i^) in order to get D > 0. This does not imply special restrictions on 
the constants involved, and, without making any particular choice for them, we can now 
write the cosmological solutions as 

a{t) = ("^'""'("^^ + "^)-"'""'^"^^ + ^)) , (3.59) 
1 B{ui sin(wit + U2) + a sin{u2t + /?)) 

(fVt) = — log —. — r — -rr^ . (3.60) 

a ^('Ui sm((X'it + ^2) — a sm(a;2^ + P)) 
As to the time behavior of the scale factor, we here limit ourselves to note that it has to 

be 

2 ' 

which soon involves a constraint on the integration constants, that is, Mi^sin^(M2) = 
a^sin^(/3). Putting U2 = (3 = in Eq. 03.591) . we soon get a(0) = 0. On the other hand, a 
choice like U2 = a = also gives the same result. Furthermore, these two possibilities both 
give a finite constant initial value for if{0). Afterwards, the scale factor evolves oscillating 
with time, describing alternate phases of closed and open universes. Whether this can be 
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the case for the universe we hve in, however, remains an open problem that needs further 
investigation. 

2. The e = — 1 value 

This situation has to be considered, now, starting from the Lagrangian in Eq. fl2.24p . 
Putting iJ^ = a'^AB then yields 

C- = u^ -v^ + 2cu\v + IrXa^iu^ - v^) + 2a^D , (3.62) 

from which we deduce the equations 

ii = -Acr^u + w^t; , V = -^a^v — bj^u . (3.63) 

It turns out that this system is of the same kind as the one in Eqs. fl3.33p . with an = 
= \XcP' , and an = —021 = 1^^, so that the solutions can be represented as in Eq. fl3.47p 
and Eq. ^Ml- 

IV. CONNECTION BETWEEN QUINTESSENCE AND INFLATION FOR EX- 
PONENTIAL POTENTIAL MODELS 

During inflation, the scalar field is supposed to be highly excited and slowly evolving 
("rolling down") to the minimum of the potential. In this phase the potential energy of 
the infiaton field dominates the energy density of the universe. After the period of inflation 
(p rapidly oscillates around the minimum of the potential and, by virtue of coupling with 
other scalar and spinor matter fields, massive and massless particles are created and the 
universe reheats, starting the standard post-Big-Bang evolution. Simple models of inflation 
use a scalar field with a potential of the form V{ip) = |m^<^^, where m is the mass of the 
infiaton field. In the standard scenario the energy of the infiaton field is then transformed 
into mass-energy of created particles. Another possibility is offered by a quartic potential 

Wi^) = a{^'-5')\ (4.1) 

where a and 5 are parameters. 

Many efforts have been made for constructing unified frameworks for inflation and 
quintessence which employ a unique scalar field to drive both stages (see, for instance, 
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Refs. 



53| 54| 55| 56|). Actually, in such scenarios the scalar field responsible of late time 



acceleration is nothing else but the remnant of the one which caused infiation at early time. 
Indeed, a successful model of quintessential infiation is subject to the constraints of both 
infiation and quintessence simultaneously. For example, the minimum of the potential must 
not have been yet reached by the scalar field (generally, this requirement is satisfied by as- 
suming the presence of a quintessential tail, that is assuming potentials with the minimum 
displaced at infinity). Moreover, it is needed that should not decay completely into a 
thermal bath of particles in order to survive until today, just to drive the late phase of ac- 
celerating expansion. As a consequence of this, the universe undergoes a period of kination 
expansion, when its energy density is dominated by the kinetic energy of (f. In this con- 
text, the standard reheating mechanism usually assumed to generate the primordial plasma 
does not work; however, the mechanism of gravitational particle production can still reheat 
the universe in the framework of quintessential infiation. Actually, even if the amount of 
radiation produced by such a mechanism is largely sub-dominant when compared with the 
energy field contribution during kination, the energy density of kination is redshifted by the 
cosmological expansion much faster than the radiation density and will start to dominate 
at some temperature. 

Generally, the constraints and requirements which are to be satisfied by quintessential 
infiation are satisfied by using a multi-branch scalar potential, where the change of the 
potential, when the field moves from the infiationary to the quintessential frame of its 
evolution, is fixed by hand 53| or is the outcome of a phase transition arranged by the 
interaction with other scalar fields (see, for instance, Ref. [57|). Recently, the possibility has 
been investigated to connect the infiationary and quintessential expansion of the universe 
within the theoretical framework of particle production, usually developed in the very early 
universe. (Even if the quantum aspect of the creation mechanism is not yet very well known, 
some classical aspects due to kinetic collision in the hot dense regions of the early universe 
have been discussed in the literature.) This suggests that the same mechanism may occur in 
ihe late universe, also leading to late time cosmic acceleration (see, for instance, Refs. jsSj 
591]). In other approaches, a large variety of quintessential infiationary potentials are derived 
from theories of non-minimally coupled gravity (see, for instance, Ref. |160i]). In our model 
the quintessential infiation is formulated in terms of a multi-branch scalar field, driving 
both the infiationary and the quintessential phases of the evolution of the universe. The 
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quintessential tail is realized through an exponential potential (hence choosing the latter as 
our working potential for quintessence); on the contrary, to describe the inflationary plateau 
we do not fix any inflationary potential, but we propose a parametrization of the inflationary 
scalar field equation of state, and implement the transition from an inflationary stage to a 
kination evolution, which is characterized by the value w^p = 1 of the equation of state, 
and corresponds just to the asymptotic in the past value for the equation of state of the 
exponential potential scalar field. 



A. The fiducial cosmological model 

In order to illustrate our paradigm, we use, as fiducial cosmological model, the one 



considered in Refs. 43| 45|]. This model is based on the simplest form of exponential 



potential of the quintessence field 

V^(y,) = ^e-Vt^, (4.2) 

and the assumption that the universe is spatially flat and filled in with dark matter and 
scalar field. We postpone to a forthcoming paper the detailed investigation of the impact 
of a kination-dominated phase generated by our class of potentials. The equations that 
determine the dynamics of this model are the Friedmann equation 

= Qm + (4.3) 

(where we use units in which SttG = c = 1, if = ^ is the Hubble constant, Qm ~ cT^ is 
the density of matter, and Q^p = \ip^ + V{ip) is the energy density of the scalar field), the 
Raychaudhury equation 

2H + ?>H' = -{^^^-V{ip)), (4.4) 

and the generalized Klein-Gordon equation describing the evolution of the scalar field 

dV 

(p + m(p + — = . (4.5) 

These equations get a simple form when instead of t as an independent variable one uses 
a(t), the scale factor. Introducing a new independent variable by m = log(l+2) = — log(^), 
where ao is the present value of the scale factor and z is the redshift, and rescaling other 
variables as = Qi = gfp- (where i = m,ip), Vq = and H = j^, where Vq is a 
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FIG. 2: The dark energy equation of state w^p parameter as a function of u. 




FIG. 3: Omega parameters as functions of u, where O^, is marked in green and fi^ in blue. 



parameter and Hq the present value of the Hubble constant, we thus obtain the following 
set of equations that contain only dimensionless variables 

,2 Qm + V 



(4.6) 



(4.7) 



with prime denoting derivative with respect to u. Such equations can be solved analytically 
but we here limit ourselves only to present the plots of some important quantities of the 
model, such as w^piu), VLmiu) and VL^{u). (See Fig. |2]and Fig. [3]for their behaviors.) 

This model of the universe is described by an exact solution of the dynamical equations. 
The arbitrary parameters that appear in the solution are determined by specifying the 
initial conditions; we in fact require that at the present time, e.g. at m = 0, one should have 
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^m{u = 0) = 0.3 and Q^{u = 0) = 0.7. The variable u is such that it decreases as time 
increases, and, at early times, the scalar field is almost constant and only recently it starts 
to increase. On the other hand, the potential of the scalar field at early times is constant 
and only recently is rapidly decreasing, while in the future it assumes a constant value again. 
The dark energy equation of state parameter in the far past is equal to —1, so that at the 
early stage of the evolution of the universe dark energy behaves as a cosmo logical constant, 
but (as is shown in Fig. |3l through the behaviors of the Q parameters) it only recently has 
started to dominate the expansion rate of the universe. 

The quintessential exponential potential admits, as said, exact solutions of the Einstein 



field equations. Actually, from them we find 42 1 



^ = -\lh^s(T^] ^ (4-8) 



3 ^\l + t\ 

which was obtained using as suitable unit of time the age of the universe, i.e. to = 1. Thus, 
at the present time (fQ = 0. At the time of reheating t is virtually zero, so that we may set 



(fin = — y 2 log (2). We now shift (p to x = if — ipin, so that the potential may be written as 

\/ = 4exp||-y|xj , (4.9) 

where the prefactor 4 is due to the particular choice of units (and, of course, is also the initial 
value for V) and represents the value of the effective cosmological constant at that time, so 
that it has the dimension of a [length]"^. In our units, the unit length is of the order of 
Hubble length, just as the current estimate for A. An effective cosmological constant near to 
4 is thus slightly greater than this and slowly evolves towards a smaller value of 2 nowadays. 

Let us now consider more realistically the inclusion of radiation, too, into such a model. 
In this case the dynamical equations as far as we know do not have analytical solutions, 
and therefore we will rely on numerical computations. Following the procedure used above, 
we use the variable u instead of time t together with the rescaled variables, so that the 
equations contain only dimensionless variables and can now be written in the form 

i^V" - {-Qr + hm + WW + ^ = , (4.11) 

2 dip 
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FIG. 4: The dark energy equation of state w^p parameter as a function of u with radiation included. 

where gr ~ is the rescaled energy density of radiation. 

We can numerically solve this system of coupled equations, specifying the initial condition 
at M = 60, for example, and assuming that v^(30), ip'{30), and H{30) have the same values 
as in the case without radiation, while ^7^(60) is just the rescaled present value of Qr- Some 
results of numerical integration are shown in Figs. H] and El It can be seen that the presence 
of radiation is slightly changing the behavior of the scalar field, its potential, the Hubble 
constant, and the w parameter of the dark energy equation of state. As expected, now, 
only the evolution of the Q parameters is different. At the initial time, again when u = 30, 
radiation dominates the expansion rate of the universe, with dark energy and matter being 
subdominant, at a redshift z of about 5000; the energy densities of matter and radiation 
become equal and, for a relatively short period, the universe becomes matter dominated 
until, at a redshift of about 1, dark energy starts dominating the expansion rate of the 
universe. With these results we then confirm the independent investigations of effects of 



radiation on the evolution of the quintessence field by Urbano and Rosenfeld ( 61[). From 
our results it indeed follows that, during the epoch of nucleosynthesis {z ~ 10^), the energy 
density of the scalar field is much smaller than the energy density of radiation. In particular, 
during such an epoch the kinetic term in the scalar field energy density vanishes, and the 
potential term is constant, so that the dark-energy term acts as an effective cosmo logical 
constant A and does not affect the process of primordial nucleosynthesis. 
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FIG. 5: Omega parameters as a function of u in the universe filled in with matter, radiation and 
scalar field, il^^ is marked in green, 0,. in red and ilm in blue. 



B. Considerations about the inflation feinaiion transition 

On the other side, as remarked above, an indispensable ingredient of quintessential in- 
flationary scenarios is the existence of an early kination-dominated (KD) era, where the 
universe is dominated by the kinetic energy of the quintessence fleld. During this era, the 
expansion rate of the universe is larger compared to its value during the usual radiation- 
domination (RD) epoch. A generic reasonable inflation-kination period can usually be de- 
scribed by two main parametrizations of the transition, without specifying the details of the 
modelization. One is a polynomial parametrization for a'^{ri) (77 being the conformal time) 



62[, and the other is a parametrization of the equation of state w, smoothly connecting the 



inflation {w = —1) and kination {w = 1) regimes. (See, for instance, Ref. 63|, where a 
hyperbolic tangent parametrization is used to study gravitational reheating in quintessential 
inflation). Moreover, the impact of a kination-dominated phase generated by a quintessen- 
tial exponential model or by quintessential power law, or also by a running kinetic inflation 



model, has been already investigated in literature (see for instance (64],[63|], [65|). Here our 
approach is quite different: actually, we are not flxing a single-branch scalar fleld potential 
which drives both the early-time inflationary and the late-time quintessence evolution of 
the universe; we instead want to show that a quintessential exponential potential tail can be 
naturally connected to the inflation through a kination dominated era, since the asymptotic 
value (in the past) of the equation of state, w^, is just = 1, which characterizes the ki- 
nation phase. We illustrate such a mechanism by using a new parametrization of the scalar 
fleld equation of state, which could also be used to study statistical properties of massive 
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non relativistic bosons arising at the first stage of relieating as a result of a quantum decay 
of a classical quantum field in inflationary cosmological model without slow rolling [66]. 
The FRW evolution of the universe in the inflationary epoch is described by the equations 



3H^ = P^, (4.12) 

(4.13) 



a I 

where we are using dimensionless units and Qq = \J^^Pi Eq^t ^ ~ 10~^^ cm. The quantities 
p<^ and P<^ are related by the equation of state P<^ = w^p^ = (7^ — 1) p^, and their evolution 
is described by the continuity equation 

p^ + 3i7(p^ + P^) = . (4.14) 

From Eqs. f|412l) . (jUS]), and flTOl) it turns out that 

7. = -?^. (4.15) 

It is straightforward to see that for 7,^ = we have an exponential behavior for a(t); if 
lip = 7o is constant, there is a power-law expansion, and we have inflation iff 70 < |, the 
more de Sitter-like, the closer 70 to zero. 

All this suggests parametrizing the dynamics of the universe during such a phase in terms 
of 7,^(t). Indeed, we have that 

m = ^ (/ l^d^ , (4-16) 
4 / . X -2 



= g (7^ - 1) ( / l^dt ] , (4.18) 



V{^{t)) = H\t) ( 1 - ^7. ) , (4.19) 



^{t) = J \l-\Hdt. (4.20) 

The above Eqs. fl4.16p and f l4.20p are parametric equations for the potential V , and can 
allow to reconstruct the scalar field potential. Thus, let us consider a scalar field which 
at the beginning describes an inflationary stage, and undergoes a phase transition into a 
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kination phase. This can be achieved by postulating a special form for (and, therefore, 
w^), and reconstructing the potential. 
Our choice for 7^ is 

where /3 and tf are parameters giving, respectively, the rate of the phase transition and 
the time at which inflation ends. On the other hand, the parameter a is related to the 
asymptotic value of w^, which in our model is w^p = 1. (For different values of a, it is 
possible to obtain a transition from vacuum into dust or radiation. Moreover, by using a 
slightly different parametrization of the equation of state, as 

= e-(^-^) + 1 " e^(-(^-t)) + 1 - ^ ' (^-22) 

) it is also possible to obtain a double transition from vacuum into kination and into radia- 
tion, as illustrated in Fig. (Q. Such a form, even if arbitrarily assigned, can easily implement 
the transition from an inflationary stage to a kination one, as we see in Figs. |6]and[71 where 
we plot the behavior of the equation of state and the reconstructed scalar field potential, 
respectively. We have to note that the potential changes rather roughly from its value during 
inflation to its final values, and that the asymptotic in the future value for the equation of 
state = 1 corresponds to the asymptotic m the past value for the quintessential exponen- 
tial potential. Thus, these behaviors can be naturally connected, as shown in Fig. |8l where 
the time-scale has been arbitrarily taken in order to show all the evolution from inflation 
into kination, toward the late-time exponential quintessential stage, and the values of the 
parameters a, /3 and tf are fixed as just optimizing the link between the early and late time 
evolutions. Finally, in Fig. [10] we compare the behavior of the scalar field potential (with 
radiation included) during the two stages: it turns out that the quintessential exponential 
potential asymptotically behaves like the inflationary one. 



V. DISCUSSION AND CONCLUSIONS 



After the historical role of scalar fields in infiationary cosmology, the quite recent recon- 
sideration of their crucial importance relies upon the fact that they can improve a dynamical 
mechanism for giving rise today to a repulsive component in cosmic energy, the dark en- 
ergy. The interest for a scalar field in the universe, otherwise usually thought of as a 
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FIG. 6: Equation of state Wip as a function of time, when the universe is dominated by kination at 
the end of inflation. The parametrization of Eq. (j4.2ip is used, with a = 2, /? = 1.32, tf = 16.9. 
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FIG. 7: Behavior of the scalar field potential as a function of time, when the universe is dominated 
by kination at the end of inflation. The parametric equations for the potential V in Eq. (j4.19p is 
used, with 7^ given in Eq. (|i:2T]) and a = 2, /3 = 1.32, tf = 16.9. 

matter-dominated one, is in fact tied to the very recent cosmological history revealed by 
astrophysical observations. 

In this paper, our attention has been generally focused on the role the exponential-like 
potential has still to play in nowadays cosmology. At first, we have discussed some models 
generalizing the simple exponential form of the potential, in order to derive general exact 
solutions. The technique used has been very simple, being it based on a given ad hoc 
change of variables. Such a procedure is generally possible either by chance or because 
there exists a sort of a method to deduce that useful transformation. In our work we 
have adopted the second procedure, i.e. the Noether Symmetry Approach to cosmology, 
borrowing from it here not only the suggestion on the kind of useful transformations to 
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FIG. 8: Evolution of the equation of state from the inflationary to the quintessential stage. Dashed 
in red the transition is plotted from the inflation into the kination when = 1; we se that such 
an asymptotic in the future value for the equation of state corresponds to the asymptotic in the 
past value for the quintessential exponential potential (solid line). 
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FIG. 9: We show the double transition from vacuum into kination and into radiation obtained by 
using the parametrization of Eq. (|i:22|) . with A = 2.1 ,B = 0.84 /3 = 8.14 ,a = -6.43 , t/ = 12.3 , 
to = 2.58. 

apply to the variables a and ip involved, but also on the natural kind of potential V{(p) 
to study. On the other hand, the results we have found in this first part of the paper are 
many and not always easily discernable. Basically, they have been mathematically derived 
but not yet appropriately discussed on physical ground. This deserves, of course, further 
investigations in a forthcoming paper, but we can however try here to sum up what has 
come out as more interesting for nowadays cosmology. 

When A = in the potential, first of all, the only case we have discussed here is for 
7^ 0,B^ 7^ 0,e = —1. The universe then displays infinitely many periods of closed 
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FIG. 10: Behaviour of the scalar field potential during the two stages: it turns out that the 
quintessential exponential potential (red dashed line) asymptotically behaves like the inflationary 
one (black dashed line). 

evolution, and one should choose which branch to investigate, anyway without succeeding 
in describing the currently observed behaviour. The situation with A > 0, on the other 
hand, depends on the sign of e, even if it presents accelerated asymptotic evolutions both 
for e = +1 and e = —1. In the latter case, such an acceleration is however obtained together 
with an oscillating behavior. This is also true for A < when e = — 1, while for e = +1 we 
again find a closed universe. 

It remains to be seen that the only physically acceptable situation (and solution) is given 
for A > and e = +1. In a forthcoming paper we are going into details of the cosmo logical 
evolution, with regard to each specific solution discussed above. In particular, we are going 
to perform the necessary confrontation of this theoretical output with the observational data 
sets. 

The second part of the paper is still connected with the exponential potential. But, 
now, we instead focus on considering and illustrating a possible quintessential inflationary 
scenario, formulated in terms of a multi-branch scalar field, driving both the inflationary 
and the quintessential phases of the evolution of the universe. The quintessential tail is 
realized through an exponential potential (choosing, thus, the latter as our working potential 
for quintessence); on the contrary, for describing the inflationary plateau we do not fix 
any inflationary potential, but we propose a parametrization of the inflationary scalar field 
equation of state, and implement the transition from an inflationary stage to a kination 
evolution, which is characterized by the value = 1 of the equation of state, which is needed 
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in any reasonable quintessential inflationary model and corresponds just to the asymptotic 
in the past value for the equation of state of the exponential potential scalar field. It turns 
out that the reconstructed potential changes rather roughly from its value during inflation 
to its final values, and the asymptotic in the future value for the equation of state = 1 
corresponds to the asymptotic in the past value for the quintessential exponential potential. 
Thus, we want to stress again that these behaviours can be naturally connected, and the 
exponential form of the scalar field potential driving the late stage of the universe could 
indeed be the asymptotic late time behaviour of the inflationary scalar field, which transits 
from inflation into kination, toward the late quintessential stage. It is worth noting that 
such a conclusion is somehow independent of the mechanism proposed for the evolution of 
the scalar field potential, being only based on a parametric description of the very early 
inflationary dynamics of the universe, and on the property of the equation of state = 1, 
which characterizes our exponential potential. In our cited forthcoming paper we are also 
going to select, among all solutions, the cases which preserve such asymptotic behaviour 
= 1. 

Of course, this mechanism for driving the transition from the inflationary evolution to- 
ward the late time accelerated expansion has to be considered as mainly exploratory, and 
some topics need more investigation. For example, the particles production, that is the 
gravitational production during the reheating or the preheating, in which particles are pro- 
duced by virtue of the variation of the classical inflaton field, needs to be investigated in a 
forthcoming paper, considering the scalar field potential given by Eq. (14.191) . (As to this, 
we have to remind that, in order to study preheating, it is necessary to couple the classical 
inflaton field ip to a massless quantum scalar field as for instance in the Lagrangian 

Ccoupied = ^d^xd'^x - l^Rx' - IgVx' , (5.1) 

where R is the scalar curvature, ^ is the gravitational couphng and 5^ is a quartic coupling 
constant between (p and x-) 

To conclude, let us note that considering radiation in the model only changes the evolution 
of the Q parameters. Radiation initially dominates on matter, while later on the energy 
densities of matter and radiation become equal; after that, for some time matter dominates 
in the universe, while dark energy starts dominating the expansion rate of the universe 
only afterwards. These results seem to confirm other investigations of effects of radiation 
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on the evolution of the quintessence field, according to which the energy density of the 
scalar field during the epoch of nucleosynthesis {z ~ 10^) is much smaller than the energy 
density of radiation. While the kinetic term in the scalar field energy density vanishes 
and the potential term becomes constant, the dark-energy term behaves like an effective 
cosmological constant, not affecting the process of primordial nucleosynthesis. 
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